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Strongly interacting quantum many-body systems are fundamentally com- 
pelling and ubiquitous in science. However, their complexity generally pre- 
vents exact solutions of their dynamics. Precisely engineered ultracold atomic 
gases are emerging as a powerful tool to unravel these challenging physical 
problems. Here we present a new laboratory for the study of many-body 
effects: strongly interacting two-level systems formed by the clock states in 
^^Sr, which are used to realize a neutral atom optical clock that performs at 
the highest level of optical-atomic coherence and with precision near the limit 
set by quantum fluctuations. Our measurements of the collective spin evolu- 
tion reveal signatures of many-body dynamics, including beyond-mean-field 
effects. We derive a many-body Hamiltonian that describes the experimen- 
tal observation of severely distorted lineshapes, atomic spin coherence decay, 
density-dependent frequency shifts, and correlated quantum spin noise. These 
investigations open the door to exploring quantum many-body effects and en- 
tanglement in quantum systems with optical energy splittings, using highly 
coherent and precisely controlled optical lattice clocks. 



1 Introduction 

Strongly correlated quantum many-body systems have become a major focus of modem science. 
Researchers are using quantum-degenerate atomic gases (1-6), ultracold polar molecules (7-9), 
and ensembles of trapped ions (10, 11) to realize novel quantum phases of matter and simulate 
complex condensed matter systems. In particular, the enhanced SU(AO symmetry in the nu- 
clear spin degrees of freedom of fermionic alkaline earth atoms may allow implementation of 
unconventional frustrated quantum magnetic models (12-15). Although degenerate gases of 
alkaline earth(-like) atoms have been achieved (16-19), reaching the extremely low levels of 
entropy required to observe these novel magnetic phases is currently not feasible in ultracold 
atom experiments. 
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Meanwhile, optical atomic clocks employing fermionic alkaline earth atoms have matured 
considerably. The most stable of these clocks now operates near the quantum noise limit (20), 
with an accuracy surpassing that of the Cs standard (27). With atom-light coherence times 
reaching several seconds, even very weak interactions (e.g., fractional energy level shifts of 
order > 1 x 10^^^) can dominate the dynamics of these systems. In fact, atomic interactions 
in optical lattice clocks were first studied in the context of density-dependent frequency shifts 
that contribute to the clocks' systematic uncertainty (22-24). When frequency shifts were first 
observed in a ^^Sr clock, they were attributed to s-wave collisions allowed by inhomogeneous 
excitation (22, 25-27), using a mean-field treatment and under the assumption that p-wave in- 
teractions were suppressed due to the ~ 1 /jK sample temperature. More recently, in an optical 
clock based on ^^^Yb atoms at ~ 10 /iK, p-wave interactions were reported to lead to two-body 
losses and density shifts (23, 28). At the same time, even at ~ 1 /^K, inelastic p-wave losses 
were observed in the ^^Sr system (29). Although the importance of many-body interactions in 
these clocks has been recognized theoretically (25-27), to date no measurements have revealed 
their many -body nature. 

In this Article, we report the first conclusive observation of beyond-mean-field many-body 
correlations in a high-density ^^Sr optical clock in a one-dimensional (ID) optical lattice. We 
perform measurements near the standard quantum limit with ensembles of > 1000 atoms. We 
show that non-classical correlations emerge as a consequence of collective elastic and inelastic 
p-wave interactions, which manifest themselves in the decay of Ramsey fringe contrast and 
as quadrature-dependent quantum noise in the effective spin degree of freedom encoded in 
the ^So and ^Pq clock states. These effects cannot be captured by a conventional mean-field 
treatment of the atomic interactions. Instead, the quantum dynamics are precisely described by 
a many-body master equation that explicitly treats quantum fluctuations in the presence of two- 
body losses. In a prior experiment (24), a strongly interacting regime (i.e., where atom-laser 
and atom-atom interactions are energetically comparable) was reached by tightly confining the 
atoms in a 2D optical lattice, at the expense of reducing the occupancy to one or two atoms 
per site. Here we probe a strongly interacting system with an average of 20 atoms per disk- 
shaped ID-lattice site, and we develop a detailed understanding of the complex many -body 
quantum dynamics. The role of s-wave collisions is suppressed by operating in the strongly 
interacting regime with highly homogeneous atom-laser coupling, making p-wave interactions, 
which operate collectively, dominant. We further show that the many-body dynamical evolution 
is well-described by a generalized infinite-range Ising model in an effective external magnetic 
field. 

The experimental observation of such quantum magnetic behavior at jiK temperatures is 
made possible because the motional degrees of freedom are effectively frozen during the clock 
interrogation. Only the internal electronic degrees of freedom (pseudo-spin) are relevant, and 
these can be initialized in a pure state. This Hamiltonian links the spin dynamics of interacting 
thermal fermions at /xK temperatures to those of two-mode Bose-Einstein condensates (BEC) 
and it has been shown both theoretically (30-33) and experimentally (34-36) to give rise to 
non-trivial many-body correlations and quantum noise-squeezed states. It is also relevant in 
trapped-ion quantum simulation experiments (10, 11, 37). The emergence of many-body cor- 
relations in optical lattice clocks subject to inelastic two-body losses opens up new scientific 
opportunities for the use of ultra-precise clocks as powerful tools for exploring strongly corre- 
lated open quantum systems (38). 
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2 Many-body /;-wave model 

We consider an optical lattice clock that employs the ^Sq — )• '^Po (henceforth \g) — )■ |e)) clock 
transition in nuclear spin-polarized ^^Sr. The lattice clock system is illustrated in Fig. [T^; it 
comprises an array of quasi-2D trap sites loaded with atoms at /iK temperatures. Fig.[l]5 depicts 
the preparation of the system including state initialization (SOM). The tight lattice confinement 
along the longitudinal direction Z freezes the dynamics and the population distribution among 
the trap sites along Z. A single site populated with N atoms is modeled as a slightly anharmonic 
2D oscillator with radial (longitudinal) frequency ur = 450 Hz (uz = 80 kHz). The motional 
degrees of freedom remain frozen to the leading order and the dynamics take place only in 
the electronic pseudo-spin (SOM). Additionally, the atoms are always prepared in the same 
state, thus Fermi statistics guarantee that no two atoms within a given trap site occupy the same 
motional state (Fig.[lll). 

We consider the aynamics within a single trap site with N atoms and label the thermally 
populated harmonic oscillator modes as nj = {nxj,nYj) with j E {1,2,... N}. Under these 
conditions, the single-site dynamics can be described by an effective spin model (SOM). For a 
pair of colliding atoms occupying two of the eigenmodes, the interactions can be expressed in 

terms of spin 1/2 operators S^'^'^ acting in the {e, g} basis as: „' ■ Sn'^ + Xnj,n'^S^^S^, + 

Cnj,nrS^/ . The subscripts rij, indicate that the spin coupling constants in principle depend 
on single-particle modes populated by the atoms. The thermal averages of the mode-dependent 
coupling constants J^, ^, , Xnj,n'., and Cn^^n' are sharply peaked about their averages J-*-, x, 
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and C with standard deviations A J-"-, Ax, and AC, respectively. However, due to the weak 
dependence of the matrix elements on the thermally populated modes, ^ A^, AC, and 
thus the large energy gap created by J-*- suppresses transitions between manifolds with different 
total collective spin S (Fig.[T^), caused by the inhomogeneities A^ and AC. Here, S{S + 1) 

is the eigenvalue of the collective operator S ■ S and S'^='-^'y^^ = ^^^^ S'^. . The dynamics of 

atoms prepared in the S = N/2 manifold is thus characterized by the collective Hamiltonian in 
the rotating frame of the laser, given by 



where constants of motion have been omitted. Here, 5 = — is the laser detuning from 
the atomic resonance and is the Rabi frequency. The terms x and C originate from elastic 
p-wave interactions and are proportional to the p-wave scattering volumes, b^g, and h^^ 
(Fig. [l]l). The matrix elements are averaged over the set of populated modes. The role of s- 
wave interactions in the absence of excitation inhomogeneity is limited to an enlargement of 
J-*-, thus enhancing the gap protection. In a 2D trap geometry, x and C are highly insensitive to 
the sample temperature T for a fixed atom number, given that the linear growth of the p-wave 
interaction with T is compensated by a corresponding 1/T decrease of the density. The term of 

order {s^^ arises when higher order corrections are used to account for virtual excitations of 

motional states (SOM). Similar effects have been measured in bosonic systems {39). Equation[T] 
is the key theoretical result of this work and describes the collective spin dynamics (Fig. [1]:). 

In addition to the elastic interaction, the effect of inelastic j9-wave two-body losses must be 
included. In contrast with ^''^Yb, where inelastic decay channels exist for both excited-excited 




(1) 
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and excited-ground collisions (23), for ^'Sr only excited-excited decay is relevant (29). The 
correct treatment of the interplay between elastic and inelastic interactions during the many- 
body dynamics is challenging due to the computational complexity of numerical techniques for 
open quantum systems. We capture the full many-body dynamics by solving a master equation 
that is numerically tractable owing to the simplification allowed by the weak {n^} -dependence 
of the p-wave interaction matrix elements. This leads to a two-body decay that preserves the 
collective nature of the model to leading order (Fig. IT]:), allowing the master equation to be 
solved efficiently for systems of up to 50 atoms (SOMJT 



3 Density shift and linear response 

Using a modified Ramsey spectroscopy sequence, we measure the density-dependent frequency 
shift of the clock transition to determine the relevant interaction parameters that characterize 
our spin Hamiltonian, Eqn. [T} The initial pulse area 9i = fiT^, chosen such that < < vr, 

controls the initial value of {S^^^}. Here, S^^^^ is the sum of over the ~ 100 relevant sites, such 

that — A'tot/S < (S^^ot) — ^tot/2, where A^tot is the total number of atoms loaded into the lattice. 

In the presence of two-body losses, {S^^^) is not constant, thus we use its time average, (5'{ot)' 

to facilitate comparison with theory. Here, we extract {S^^^) from independent measurements 
periodically inserted into the clock sequence. The duration of the dark time, Xdark, is fixed at 
80 ms and the final pulse area is set to 7r/2. We measure the shift by modulating the density by 
a factor of ~ 2 (Fig.JZl). 
A simple mean- fie 



d analysis of Eqn. [T] (neglecting cubic terms and losses) treats interactions 

as an effective magnetic field along Z, B{N) = NC + 2x{S'), where (S^) = -{N/2) cos^i. 
The mean-field density-dependent frequency shift Au^N) = B{N)/ (27r) scales linearly with 
the excitation fraction and agrees with experimental observations (Fig. [2|. Additionally, we 
fit an exact solution of Eqn.|l]to the data. Both fits are shown in Fig.lzj To compare with 
the experiment, we always perform an average over the atom number distribution across the 
lattice sites. From this measurement we extract x = 27r x 0.20(4) Hz and C = — 0.3x- The 
linear behavior of the density-dependent frequency shift with Ramsey spectroscopy was also 
measured in the ^^^Yb lattice clock (23, 28). 

As a complementary measurement, we further probe the linear response dynamics of the 
system spectroscopically with a single weak pulse while varying 5. The pulse duration is Tr = 
500 ms, and its area is 0,Tr 0.27r. At this low excitation fraction, two-body losses can 
be ignored, and the mean-field analysis yields a density shift of the clock resonance given by 
Az/LRR ~ (C — x)N, which agrees with the mean-field expression of the Ramsey frequency 
shift at small 6i. In the experiment, we observe a shift of 2.7 Hz for a modulation between 
-^High = 3 X 10'^ and A^low = 1 x 10'^ (Fig. |2| inset). We extract x from a simple mean- 
field model and find agreement at the 15% level, consistent with the atomic number distribution 
uncertainty. 



4 Rabi spectroscopy 

Having extracted the interaction parameters as described in Section |3} we turn our attention to 
regimes where strong interaction effects emerge by probing the system via Rabi spectroscopy. 
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The atom-laser interaction energy, h^l, is reduced to a level comparable to the many-body in- 
teraction energy. The pulse duration, Tr, becomes correspondingly long to allow for high 
spectroscopic resolution of the many -body excitation spectra. Typically we set ^ITr = tt. In 
terms of the spin Hamiltonian, the evolution now takes place in the presence of a transverse 
magnetic field-like term that does not commute with the Ising interactions, thus giving rise to 
non-trivial quantum behavior. In a direct analogy to a double-well potential for BEC, and 
S play the role of tunneling and bias between the left and right wells, respectively; and their 
competition with non-linear interactions has been shown to lead to novel quantum dynamics 
including macroscopic quantum self-trapping (40,41). 

At the highest achievable densities in our experiment, we observe dramatic deviations in 
the Rabi lineshapes from ideal, single -particle lineshapes. By varying either the atom number 
or n (with Tr = vr/fi), we change the ratio of the many-body interaction energy to hVl. The 
resulting lineshapes are summarized in Fig.|3| Our single-scan frequency resolution is 500 mHz, 
and we are thus capable of resolving sub-Hz mteraction features spectroscopically (Fig. [3^). At 
low atom number (A^tot — 1 x 10^) and with Tr = 200 ms, we observe a nearly ideal Rabi 
lineshape that becomes severely distorted with increasing density (Fig. Similarly, as we 
increase Tr, keeping i^T/j = tt, we begin to see the onset of an interactionblockade mechanism 
(macroscopic self-trapping in the BEC context (40)) in concert with the distortion (Fig.[3j:). At 
the largest Tr = 750 ms, we excite only 20% of the atoms into |e) and the line is approximately 
five times broader than the non-interacting lineshape. 

Using the interaction parameters extracted via the Ramsey density shift and the linear re- 
sponse measurements, we are able to theoretically reproduce all the observed lineshapes (the- 
oretical curves shown in Fig. |3b-c) using entirely the mean-field treatment of Eqn. [T]in the 
presence of two-body losses. Here, a full many -body treatment of the master equation is in 
agreement with the mean-field predictions. 



5 Quantum fluctuations and beyond mean-field physics 

While the density-dependent frequency shift and lineshape measurements provide verification 
for the mean-field limit of the many-body model and are useful for extracting the interaction 
parameters in a simple and efficient way, they do not represent a sensitive probe of the full 
many-body dynamics. Ramsey fringe contrast, on the other hand, can undergo a periodic series 
of collapses and revivals, which will reflect the quantized structure of the many-body spectrum 
arising from the granularity of the interacting atoms. Specifically, for a fixed and neglecting 

both the losses and the (S^)^ term, the amplitude of the Ramsey fringe evolves as Z^^^ with 
= 1 — sin^(6'i) sin^ (x^dark)- At short times (i.e., V^X^dark ^ 1)^ the contrast is approxi- 
mately 1 — (A^ — 1) sin^(^^i)r|jjj.]^/2. On the contrary, the mean-field model at the single-site 
level (with fixed A^) predicts no decay of the Ramsey fringe contrast, as a magnetic field B{N) 
can only lead to a pure precession of the collective Bloch vector. By taking the average over 
atom distributions among lattice sites and properly treating two-body loss during the Ramsey 
dark time, the mean-field model does show a decay of the contrast. However, this decay is 
associated mainly with dephasing arising from different precession rates exhibited by sites with 
different A^. This effect is most relevant for low excitation fractions where B(N) is large and 
the inelastic loss is suppressed. 

We implement a Ramsey sequence to experimentally measure the fringe contrast. The Ram- 
sey pulse durations are < 6 ms, satisfying Cl ^ Nx, to suppress interaction effects during the 
pulses. We apply the final tc/2 readout pulse with a variable relative optical phase of 0° - 360°. 
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We record the fraction of excited atoms as a function of the readout phase, which constitutes 
the Ramsey fringe. The Ramsey fringe contrast is extracted in a manner that is insensitive to 
the frequency noise of the ultrastable clock laser (SOM). 

We explore three distinct experimental conditions to rule out single-particle decoherence 
mechanisms and thoroughly test the model. The first condition represents the typical operating 
parameters of the lattice clock, with A^tot = 4 x 10'^ and vz = 80 kHz. In the second case, 
we reduce the lattice intensity such that = 65 kHz, which results in a reduction of the 
density by a factor of ~1.8. Finally, we maintain uz = 80 kHz but reduce the atom number 
to A'tot = 1 X 10^. Under all conditions, the full many-body density matrix model reproduces 

the experimental observations well (Figs.|4^,|4|:, and 4^). The inclusion of the {S^Y correction 
improves the theory-experiment agreement, especially for pulse areas 6i > 71-/2 and for the 
high-density conditions (SOM). We also observe a striking breakdown of the mean-field model 
for 9i > n/2 where many-body corrections are expected to be dominant. 

The frequency shift, lineshape, and Ramsey fringe contrast are quantities that all depend 

on the first-order expectation values of the spin operators (S*^'^'^). We now turn our attention 
to the distribution of quantum noise, which depends on the second-order moments of the spin 

operators, e.g., { (^sA ) - (S^)^, {S'^S' + S'S'-") - 2(5^') etc. Given that the form of the 

Hamiltonian in Eqn.[T]is known to produce squeezed and entangled states (31), the distribution 
of the quantum noisebecomes a compelling measurement to probe physics beyond the mean 
field. However, until this work, optical clocks have lacked the requisite laser precision to probe 
the quantum noise distribution. Recent advances in ultrastable lasers have permitted access to 
the quantum regime (20). 

To minimize single particle decoherence effects and dephasing due to the distribution of site 
occupancies, we add a spin-echo pulse to the Ramsey sequence. As a result, the sensitivity to 
low-frequency laser noise is reduced at the expense of increased sensitivity to high-frequency 
laser noise. With atoms initialized in \g), we follow the pulse sequence shown in Fig.lSlto ma- 
nipulate and measure the spin noise of the many-body state. For each value of the finalrotation 
angle, representing a specific quadrature in which we measure the spin noise, we repeatedly 

record {S^^^)/Ntot via measurements of the final atomic excitation fraction after the Ramsey 

sequence. From the data, we determine = {(^S^^^^ )/^tot ~ ('S'tot)^/^tot analyzing the 

pair variance for successive measurement of (S^^^) /A^tot- We note that the quantum limit of cr^ 
is important for defining the ultimate stability of lattice clocks (20). For an ideal coherent spin 
state of the entire ensemble, the standard quantum limit of is given by cr^qj = pi^ — p) /^tot> 
where p is the probability of finding an atom in the excited state, and can be estimated as 

p={S^)/Ntot + l/2. 

We perform measurements under three conditions in order to probe the time evolution of 
the quantum noise distribution: Tdark = {20 ms, 40 ms, 60 ms}. For these dark times, we used 
TT -pulse times of Tr = {10 ms, 20 ms, 20 ms}, respectively. We find that utilizing long vr pulses 
reduces the sensitivity to spurious high-frequency components of laser noise. As summarized 
in Fig. [sj for Ntot = 1 x 10'^, the quantum noise contribution to the spin noise is comparable 
to that of the laser noise. However, with Ntot = 4 x 10^, the laser noise is responsible for a 
larger fraction of the noise in repeated measurements of (S^). There are qualitative differences 
between the low and high atom number cases; for example, for Ni^i = 4x10^ with Tdark = 20 ms 
and 40 ms, we observe a phase shift for the minimum of the spin noise. We have verified that 
in the limit of small Tdark, the phase shift is no longer prevalent. To compare the predictions 
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of the full many-body master equation with the experiment, we add the effect of laser noise 
in quadrature with the calculated spin quantum noise. We additionally find it necessary to 
fully treat the effects of interactions during the laser pulses. As shown in Fig. |5| the theory 
predicts a direction and magnitude of the phase shift of the noise minimum that agree with the 
experimental observations. The theory also predicts significant spin noise for rotations near 
±90° due to an effect analogous to anti-squeezing. Despite the high sensitivity to laser noise at 
this rotation quadrature, the measurements of the total noise are consistent with theory. 



6 Conclusion and outlook 

The use of an intricate competition between coherent interactions and dissipation as a means to 
enhance and engineer many-body correlations is an emerging theme in modern quantum science 
(42^4). We expect that the capabilities demonstrated here will develop into important tools to 
explore the rich interplay between many-body effects and dissipation in alkaline-earth atoms 
confined in optical lattices, as well as in other systems such as polar molecules (7-9, 45, 46). 
Furthermore, as a matter of practical importance, these techniques could be used to generate 
interaction-induced spin squeezing in an optical lattice clock {38). Although the investigation 
described here is restricted to nuclear- spin-polarized gases, exploration of similar many-body 
effects in a clock making use of additional nuclear spin degrees of freedom is a promising first 
step towards the investigation of SU(A^) orbital magnetic models in ultracold atoms. 
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Figure 1: Schematic diagram of the interacting many-body system, (a) Approximately 100 
lattice sites are significantly occupied during the experiment. The average lattice occupancy is 
20 atoms, (b) At t = 0^, the particles are initialized in the ground state. After a short pulse, 
the atoms can be prepared in a product state of arbitrary single particle states, {a\g) + I3\e))'^^ 
at t = 0+. (c) The product state experiences the all-to-all many -body interaction for t > 0. (d) 
Interactions between particles are parametrized by the mode- and spin-dependent interaction 
parameters, Vgg oc 6^^, Veg oc 6^^, and Vee oc (see SOM for definition), (e) The many- 
body Hamiltonian has eigenstates comprised of maximally symmetric superpositions (Dicke 
states, for which S = N/2) of electronic ground and excited state superpositions, depicted as 
purple shells. Because of the slight inhomogeneities in the coupling strengths, the maximally 
symmetric manifold is coupled to the next lowest manifold with S = N/2 — 1, depicted as a 
nested blue shell. However, this coupling is prevented by an energy gap present in the many 
body Hamiltonian as a result of the J-^S ■ S term, which is depicted by an offset of the two 
manifolds in this figure. Two-body inelastic losses connect maximally symmetric manifolds of 
S S — 1, and thus are not a strong decoherence mechanism. 
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Figure 2: Density shift in Ramsey spectroscopy fit with the full many-body solution. Due 

to the perturbative nature of this measurement, the mean-field approximation to the many-body 
theory (dashed line) agrees well with the data. The exact many-body solution in the absence 
of losses (solid curve) agrees best with the data only for lower values of total average spin 

(S'tot) /Ntot due to the nonlinear (S*^)^ term in the Hamiltonian. The zero crossing occurs at an 
average excitation fraction, given by (Sfot) /^tot + 1/2, of = 0.64(1). From the zero crossing 
and the measured slope we extract x ^iid C. At (5't^ot)/^tot — ±1/2, the dephasing due to 
population differences in sites is maximized, resulting in maximal Ramsey fringe decoherence 
(shown in Fig. |4]). Near the zero-crossing, the Ramsey decoherence is dominated by intra-site 
effects due to the nonlinear term in the many-body Hamiltonian. Inset: linear response of the 
Rabi spectroscopic line shapes under weak excitation (pulse area 6i = 0.2n). The shift of the 
line center provides an independent verification of the model in the linear response regime. 



9 




-6 -4 -2 2 4 "-^ 50 100 150 200 250 300 



Detuning (Hz) Detuning (Hz) 

Figure 3: Lineshape as a function of experimental parameters. Each curve is a composite 
of multiple scans that have been centered atop one another. These data are subsequently binned 
with a binsize of 1 Hz. The curves are offset in both the vertical and horizontal directions for 
visual clarity, (a) Lineshape obtained at an operating density of O.lpo with a 2 s clock probe, 
where po is the density obtained for 5 x 10'^ atoms, and is given by po = 5(2) x 10^^ cm^^. 
The line shape shows a clear distortion due to the density-dependent interaction at the sub-Hz 
energy scale. (Inset) Extremely low density scan ( p = 5 x lO^^po) with a 3 s clock probe, 
demonstrating the ultimate frequency resolution of the system, (b) Lineshape as a function 
of density (normalized by po)' with a 7r-pulse time of 200 ms. (c) Lineshape for p ~ po as 
a function of probe time. In both (b) and (c), theoretical curves, obtained with the mean-field 
treatment including loss, are shown as dashed black lines and agree well with the measurements. 
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Figure 4: Ramsey fringe contrast decay vs initial pulse area. The pulse area values for 
the data and corresponding theory are given in the legends of the plots. Error bars represent 
the statistical error of each contrast measurement, and thus do not account for systematic drifts 
that occur over the course of the experiment. The solid lines (top panels) are the many-body 
calculations, while the dashed lines (bottom panels) are using the mean-field approximation of 
the theory. The many-body model and the mean-field approximation agree in the limit of small 
initial pulse area {i.e., Bloch vector polar angle), but disagree for pulse areas > 7r/2. This is 
an important confirmation of the dominance of many-body effects in this parameter regime, (a, 
b) uz = 80 kHz, ur = 450 Hz, and iVtot = 4000; (c, d) uz = 65 kHz, ur = 370 Hz, and 
iVtot = 4000; and (e, f)uz = 80 kHz, ur = 450 Hz, and A^^tot = 1000. 
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Figure 5: Spin noise vs. quadrature, (a) An initial pulse prepares a coherent state along y, 
which then evolves for rdaik/2. An echo pulse then rotates the many-body state 180° about — x. 
After and additional evolution time of rdark/2, a final pulse rotates the state about —y and the 
spin noise is measured. The many-body state depicted here represents the spin evolution of a 
20 atom ensemble in a single trap site with r^ark = 40 ms. To remove spurious effects due to 
slow drifts in atom number, the data is processed as detailed in the SOM, which removes the 
potential bias. For plots (b-g), the dashed line is the pure laser noise extracted from a fit to the 
data. The solid line is the laser noise plus the full many-body prediction of the spin noise. This 
full theory is simultaneously fit to both the low and high atom number curves to extract the laser 
noise for a given dark time. Vertical arrows indicate significantly phase-shifted minima in the 
experimentally measured spin noise, consistent with the predictions of the many-body theory. 
Plots (b, c): Spin noise for A^tot = 1x10^ and Atot = 4x10^, respectively, with an interaction 
time of Tdark = 20 ms and a vr -pulse time of 10 ms. A clear phase shift is observed near the 
origin for the high atom number case, (d, e): Spin noise for Atot = 1 x 10'^ and Atot = 4 x 10^, 
respectively, with Tdark = 40 ms and a vr -pulse time of 20 ms. In plot (d), a clear phase shift 
is also observed near the origin, (f, g): Spin l^ise for Atot = 1 x 10'^ and A^ot = 4 x 10^, 
respectively, with Tdark = 60 ms and a vr -pulse time of 20 ms. 



Supplementary material 



Experimental procedures 
Atom preparation 

We load fermonic ^"^Sr into a one-dimensional vertically oriented optical lattice via two- stage 
laser cooling, directly producing lattice-trapped samples with temperatures of 3-5 fiK. The op- 
tical lattice is maintained near the magic wavelength (47, 48) for the ^'^Sr ^Sq ^Po clock 
transition. After further cooling and optical pumping into ^Sq |-F = 9/2, mp = +9/2), we 
obtain 5 x lO'^ nuclear spin polarized atoms concentrated mainly in a central region of approx- 
imately 100 disk-shaped lattice sites with typical trap frequencies as described in the text. We 
determine the longitudinal temperature, Tz, via sideband spectroscopy and the radial tempera- 
ture, Tr, via Doppler spectroscopy (24,49). With typical sample temperatures of = 1—2 fxK, 
Tr = 2 — 4 ^K, we find that for our trapping conditions the average density is 5(2) x 10^^ cm~^ 
for a sample consisting of 5 x 10^ atoms. With the sample prepared as outlined above, we per- 
form high-resolution spectroscopy on the ^So l-F = 9/2, mp — +9/2) ^Pq \F — 9/2, mp — 
+9/2) clock transition. 

Contrast measurement 

As described in the text, we measure Ramsey fringe contrast as a function of initial pulse area. 
For a first pulse of area 9i, we allow the system to evolve for time r. We then apply a final pulse 
of area 7r/2 and measure the resulting excitation fraction as a function of the optical phase of 
the second pulse relative to the first pulse. For r > 100 ms, there is a significant additional 
random phase added due to the frequency fluctuations of the ultrastable clock laser. A given 
excitation fraction (pi) measurement will yield Pi = Csin^ where C is the contrast and 

A0j is the ith realization of the both deterministically and randomly varied phase. By analyzing 
Var (p) = and assuming a uniform distribution of A0i, we obtain the contrast in a way 
that is insensitive to the laser noise. 

Data analysis for spin noise measurement 

We perform quadrature-dependent spin noise measurements as described in the text at a given 
target atom number A^tot = 1 x 10^ or A^tot = 4 x 10'^ During the course of these measure- 
ments, we typically observe slow, systematic fluctuations of the atom number on the order of 
±10% as we operate the experiment and measure spin noise over the course of ~10 hours. In 
most instances, these fluctuations are negligible due to the normalization techniques we employ. 
However, the atomic spin noise depends directly on the atom number, and a slowly varying atom 
number could result in unintended systematic biases. Specifically, spin noise for the coherent 
spin state typically considered in optical clocks scales as 1/ V^tot- Thus, the deviations in atom 
number can cause variations on the order of ±5% in the measured spin noise. Ideally, these 
fluctuations are randomly distributed and should not result in interpretation as a false-positive 
for non-trivial spin-noise correlations. In the unlikely possibility that these fluctuations were 
correlated with a specific measurement quadrature, they could cause a spurious phase shift in 
the spin noise minimum. We thus take care to analyze the data in a way that is immune to 
this potential bias. In this supplement, we describe our process for removing the variability in 
the spin noise data due to a slowly fluctuating total atom number. In this way, we verify that 
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the phase shifts measured in Fig. 5 are not manifestations of the "trivial" case, where the spin 
noise is described by a coherent spin state. Rather, we verify that the phase shift of the spin 
noise minimum is a direct result of the many-body interaction as described in the Article and 
Methods. 

A given measurement of {S^^^)/Ntot is accomplished by independently measuring N(,(g), 
the number of atoms in the excited (ground) state after a single Ramsey experimental sequence, 
using standard electron shelving techniques. We determine its ith value, (>S'tot)i/^tot' t'Y mea- 
suring the ith value of -/Ve(ff) (which we denote as ^l^^-^) and obtain 

{SL)i/NL = -J^, - 1/2. (2) 

From the jth set of measurements of 0^^^), denoted {{S^ot)i^ ■ ■ ■ > {^tot)h ■ ■ ■ {Stot)nj}j, we 
estimate = {(S^j'^)/Nl^ - {S^J'^/Nl^ using a pair variance, such that 



^E(('^tot)m-(^tot).)'- (3) 



For white noise, the pair variance is a good estimator for the standard deviation (50), while 
remaining insensitive to noise processes that only manifest themselves on long time scales. 
The number of measurements in a set, nj, was typically nj ~ 80. For a given measurement 
quadrature, we average the results of many such measurement sets to produce one experimental 
data point (i.e., a data point in Fig. 5). 

In order to maintain insensitivity to slow fluctuations in atom number between sets j and /, 
we consider the standard expression for quantum noise for the case of a coherent spin state, crgqi, 
which is expected in the absence of many-body interactions. The explicit goal is to remove any 
mechanism by which the trivial case — where the spin noise is described by crgqi — can mimic the 
many-body effect we predict from the theory. We calculate the jth value of (Tsqi as 

(<i)'=P.(l-Pi)/<t, (4) 

where pj = Mean [{N^/ {N^ + iVj) , . . . , N^/ (Ni + N^),..., Ne' / {K' + Ng')}j] . We 
additionally consider a technical noise term, which represents the effect of intrinsic technical 
detection noise, given by Asj. This noise is characterized by a separate measurement. The 
detection noise accounts for 10% of the observed noise at typical low atom numbers, while 
at high atom number it is only ~ 1% of the observed noise, and is therefore negligible. It 
is quadrature-independent in all cases. From the cr|, we subtract the atom-number-dependent 

«i) such that 

^' = '^-«i)'-A4 (5) 

Here, represents only the effects of non-trivial spin noise and laser noise. 

The many-body theory for a given measurement condition is calculated at fixed atom num- 
ber. To facilitate comparison with the many body theory, we add a noise term back to that 

corresponds to cr^qj for the mean atom number over the entire data set, cr^qj. We emphasize that 
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(7gq, is a constant number, with no quadrature dependence. The many-body theoretical predic- 
tion is calculated based upon the same mean atom number used to calculate 0-^^,. Ultimately, the 
net effect of this process is to remove the variability due to slow fluctuations in atom number, 
but to retain the part of the noise that departs from cxsqi due to both laser noise and many-body 
effects. As discussed in the text, we observe a phase shift of the minimum of the phase noise 
that is consistent with the many-body theory and indicative of correlated spin noise of the atom 
ensemble. 



Derivation of the spin Hamiltonian 

The many-body Hamiltonian describing a nuclear spin-polarized ensemble of fermionic atoms 
with two accessible electronic states, g and e, which experience the same external potential 
V(.xt{^, can be expressed as 



zm I m ' ^ 



m 



a,/3 L J L 



(6) 



Here ^a(R) is a fermionic field operator at position R for atoms with mass m in electronic 
state a = g,e. We have included only s-wave and p-wave channels, an assumption valid at 
/iK temperatures. Since polarized fermions are in a symmetric nuclear-spin state, their s-wave 
interactions are characterized by only one scattering length a~g, describing collisions between 

two atoms in the antisymmetric electronic state, :^(|fi'e) — \^9))- The p-wave interactions 

can have three different scattering volumes 6^^, fe^g, and h^^, associated with the three possible 

electronic symmetric states ( \gg), |ee), and + \Gg)), respectively. Note that here we 

have assumed no external laser field and have ignored the optical energy splitting since the 
number of atoms in the excited state remains fixed in the absence of driving terms. 

In the experiment, V^xtiR) is a deep ID lattice along Z, which creates an array of two- 
dimensional discs and induces a weak harmonic radial (transverse) confinement with an angular 
frequency ujr = 2txur. The lattice confines the atoms to the lowest axial vibrational mode. We 

expand the field operator in a harmonic oscillator basis, ^q(R) = <Po{Z) X]n^"n</'nx(^)*/'"y(^)' 
where 0f and 0„ are, respectively, the longitudinal and the transverse harmonic oscillator eigen- 

modes and cj^^ creates a fermion in mode n = (n^, riy) and electronic state a. In this basis, H 
can be rewritten as (23, 24) 



h 



a,n 



a,/3,n,n',n",n" 



-f, = 



"ho 



(7) 



Sa,i3 is a Kronecker delta function. Here, af^ = ^/h/ {mujji) is the radial harmonic oscil- 
lator length, and 5'nn'n"n"' and Pnn'n"n"' characterize s- and p-wave matrix elements respec- 
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^nn'n"'n" 



lively which depend on the harmonic oscillator modes and satisfy S'nn'n"n"' = S„ 

5'n'nn"n"' = 'S'n'nn"'n" and Pnn'n"n"' = " -Pnn'n"'n" = " -Pn'nn"n"' = -Pn'nn"'n"- -E'n are singlc- 

particle energies in the trap. 

Under typical operating conditions, uji ~ 450 Hz, the interaction energy is about two orders 
of magnitude weaker than the single -particle energy. Thus, at the leading order, only collision 
events that conserve the total single-particle energy need to be considered. For a slightly anhar- 
monic spectrum (the anharmonicity is weak but can be at the level of the interaction energy), 
with only one exceptiorQ those processes conserve the number of particles per mode. In this 
case, for an initial state with at most one atom per mode (| (7) -polarized state), it is possible to 

reduce H toa spin- 1/2 model. 

N Q 

H'/h = Y}Jtr.,XS^. ■ Sn,) + Xn„n,,^n,^;., + + ^^0]- (8) 

Here S^^ = \ p ^anj^afsc^nj, with cr^'^'^ Pauli matrices. Constant terms have been dropped. 

V^9 _ Jjeg yee , ygg _ 2\/^9 (V'' - V^^ ' 

nw.n / o ' Xn,,n./ „ ; '-^n,,n./ „ 



The quantities 

encapsulate the temperature dependence of the interactions. A further simplification of Eqn.[8] 
can be made thanks to the following considerations: (i) The atoms start in the totally symmetric 
manifold with S = N/2 since at time t = all the atoms are polarized in the g state; (ii) 
P-wave interactions are suppressed by the centrifugal barrier, and thus, at /xK temperatures, U 
should in general dominate over V (although it must be said that the actual values of the s and 
p scattering parameters are not known); (iii) The p-wave matrix elements, Pnj,n-,, are functions 
which weakly dependent on rij , n^/ . All of these considerations generate an energy gap that 
prevents transitions between the S = N/2 and S = N/2 — 1 sectors. Consequently, to a very 
good approximation the dynamics can be projected into the S = N/2 manifold with an effective 
Hamiltonian given by Eqn. 1 (which also contains single -particle terms arising from the optical 

driving field), x = ^%^(n-i)^' ^^'^ ^ ^ ^^NiN-i)"^' mode-averaged quantities. We have 
omitted the term S ■ S, which is a constant of motion. 

Decoherence and two-body losses - Exact treatment 

The Hamiltonain formulation described above is valid only for a closed system. To account for 
losses due to inelastic e-e collisions, one needs to use instead a master equation: 

hj^p = -i[H,p] + Cp. (11) 

'The only exception are the processes (nx^ , ny^ ) and (nxa i "'y2 ) ~^ ('^^2 1 ) and {nx-^ , n-Ya )• 



16 



Here p is the reduced density matrix operator of the many-body system. H is the Hamiltonian 
given by Eqn. 1 , and £ is a Liouvillian that accounts for inelastic processes. Considering p-wave 
e-e losses and under the same assumption of frozen motional degrees of freedom, C is given by 

Here the jump operators are An, n , = Cen,Cen , and Fn, n , = 7'^'^-Pn, n ,• The expression for 7 

■'3 •' 3 ■'3 ■'J 

is identical to v^'^ up to the replacement of the p-wave elastic scattering volume by the inelastic 
one. 

An important note is that there is no coherence between sectors of different atom numbers 
and the master equation can be solved in a "block-diagonal" way (see Fig. 1). Specifically, if 
the system starts with N particles, we need to solve a series of differential equations for each 
of the subspaces with cascading atom numbers: first for N particles, next for — 2 particles, 
then for — 4 particles, etc. There are (^) = n\{N-n)\ different sectors with N — n particles, 
n = 0,2 ... N. For mode independent interaction parameters ^nj,n^, T. In each of the 
J\f = N — n sectors, the dynamics is restricted to the collective Dicke states \S = M/2, Mj^) = 
\Mj^). Moreover, each of the (^) sectors behaves identically. 

Let pj^ be the density matrix for a single sector of J\f particles. Furthermore, let us assume 
that particles in are numbered from 1 to A/" in such a way that atoms J\f + 2 and J\f + 1 are 
the ones that decay as one goes from pAr+2 to p_^f. The resulting equations are 



i<j ^ ' 



for < A/" < A^. In terms of spin operators. 



+ + (14) 



(m ^ 2) (m' ^ 2) 

(MAr|Ce,Ar+2CeX+lPA^+2Ce,Ar+2^eX+ll^A^) = (^A^+2 + '^\PAf+2\M^^^ + 2) ^^p^ w y '^^^^ 

For < 50, the above equations can be efficiently solved numerically. |^ 

Decoherence and two-body losses - Mean-field-treatment 

A simple and illuminating way to perform a mean-field treatment is to use the Schwinger bo- 
son representation that maps spin operators to two-mode bosons (51). It represents the spin 
operators as 

2S' = ^1^1 - ^J^o S+ = ^1^0 S- = (16) 

^ At tfie highest operating densities a cut-off of 50 atoms per lattice site is sufficient. 
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with a bosonic annihilation operator of mode j = 0, 1. In terms of Schwinger bosons, the 
collective spin model maps exactly to the Hamiltonian that describes a double-well bosonic 
system (by identifying the and 1 states with the left and right wells) with a tunneling matrix 
element ^l, a bias between the wells 5, on-site interactions V^^, V^^ in the left and right wells, 
respectively, and nearest-neighbor interactions V^^. The mean-field approximation, which gives 

rise to the Gross-Pitaevskii equation, replaces the field operators o with complex numbers, 

^1,0 ^0,1- To deal with the lossy dynamics, it is necessary to express the equations of motion 
in terms of the density matrix, pij = The resulting equations are 



d_ 

Ft 
d_ 

Ft 
d_ 

Ft 



Poo = -i^iPio- Poi), 
Pu = i^ipio- Poi) - 
Pio = ^ttIPii - Poo) + « 



rp?i, 



-6 + C{pii + Poo) + x(Pii - Poo) + i 



rp 
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(17) 
(18) 
Pio, (19) 



where poi 
1. 



p*Q . The mean-field treatment is expected to be valid in the regime v sin [6i] xt ^ 



Higher Order Corrections 

The spin Hamiltonian neglects collision processes that do not preserve single -particle energy. 
However, collisions populating off-resonant states can still take place virtually and will intro- 
duce corrections to the spin model. To account for those, we split the Hilbert space into the 
resonant, S, and non-resonant, T, manifolds respectively, spanned by the states 

N 



= |a„„a„„...,a„^), E*"*^^^^^^., (20) 

i=i 

N 

\^\) = ^ J]i5k, (21) 

j=i 

These states are written in the occupation basis and a E {g, e}. In lak^ ), the same mode kj can 
also be occupied by both an e and a g atom simultaneously. 

By projecting the many-body Hamiltonian on S using the corresponding projection operator 

one can obtain an effective Hamiltonian: 

= V^HV^ = + H^^ . (22) 

Here, is the spin model given by Eqn. [s] H^^ can be obtained via the second order pertur- 
bation theory as follows: 

($5 1^52 1$^ N _ _ V ^ "g^^ "fe"' (23) 

ar k 
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When H^'^ is written in terms of spin operators and is projected into the collective Dicke 
manifold, it gives rise to terms proportional to S^,{S^Y and (5*^)^. We note that terms propor- 
tional to {S^Y vanish for p-wave interactions. Among those corrections, the first two can be 
absorbed into the spin Hamiltonian given by Eqn. 1 and only the third one gives rise to addi- 
tional corrections. The corrections in the Ramsey contrast decay arising from the cubic term 
are shown in Fig. S 1 . 
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Figure 6: Ramsey fringe contrast decay vs initial pulse area. The solid lines are the many- 
body fits using terms of order {s^^ , and the dashed lines represent fits without these terms. 
These data are the same as in Fig. 4a of the main text, and were taken with vz = 80 kHz and 
iVtot = 4 X 10^. Pulse areas are shown in the legend. 
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